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Abstract 

We estimate the numbers and mass spectra of observed lepton and kaon pairs 
produced from (j> meson decays in the central rapidity region of an Au+Au 
collision at lab energy 11.6 GeV/nucleon. The following effects are consid- 
ered: possible mass shifts, thermal broadening due to collisions with hadronic 
resonances, and superheating of the resonance gas. Changes in the dilepton 
mass spectrum may be seen, but changes in the dikaon spectrum are too small 
to be detectable. 
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I. INTRODUCTION 



There is currently a great deal of interest in finite-temperature properties of hadronic 
resonances Recently, preliminary results for the (f) meson mass spectrum in Si+Au 

collisions at lab energy 11.6 GeV/nucleon were obtained at Brookhaven's Alternating Gra- 
dient Synchrotron (AGS), by reconstructing the <fi mesons from final-state kaon pairs ||. 
No change was observed from the vacuum width of the 0, although a possible small mass 
shift was observed for the most central events. These negative preliminary results could 
dampen the enthusiasm of other groups to study the <fi peak in the dilepton channel. We 
show here that no observable change is expected for the <fi peak in the dikaon spectrum, 
but that effects may be visible in the dilepton spectrum. Thus, there is still good reason to 
study the dilepton peak at the AGS, in spite of the fact that dikaon results are negative. 



II. EQUATION OF STATE 

The behavior of the hot matter is somewhat more complicated in events at AGS energies 
than in ultra-relativistic events, due to the large baryon densities and the lack of strangeness 
equilibration. The equation of state used here is almost the same as those used in Refs. |9|-[TT 



We describe the system using the temperature T, the baryon chemical potential \1b, and 
the strangeness and antistrangeness chemical potentials, respectively [is and fig, using the 
high energy conventions c = H = ks = 1- 

The matter has nearly isospin zero, so the up and down quark chemical potentials, 
respectively \i u and fid, are fi u = fid = A*b/3. The u and d quarks are approximately in 
chemical equilibrium with their antiquarks, so the u and d chemical potentials are /% = 
fjq = — /is/3. Finally, the net strangeness is also zero, since strong interactions conserve 
strangeness and the time scales are too short to allow weak interactions to be significant. 
The s and s chemical potentials are respectively fi s = Hs + Hb/S and fig = fi-g — /-ts/3, so that 
H~s = A i 5 + 2/i J B/3 [fTIfl . It is unlikely that s and s are in chemical equilibrium with each other, 
because VR S <C n s dV/dt, where V is the volume of the hot matter, R s is the production 
rate of s quarks per unit volume, n s is the density of s quarks, and t is time. However, we 
assume for simplicity that the system is in chemical equilibrium, so that fis = —/is/3; this 
assumption will be relaxed in later works. The mean number of s quarks per event is much 



greater than unity, so the possible ss pair chemical potential [12| is omitted. 

We model the QGP as a collection of free quarks and gluons, with a constant (bag) 
energy density B ~ (220MeV) 4 . The low-temperature phase is treated as a resonance gas 
(RG), using all confirmed strongly-interacting particles with known quantum numbers and 



without c or b content JDJ. The chemical potential for resonance i is 

tk = + + A« - A?> - A® - Af) /V3 + A?W + A?W (1) 

where is the number of constituent quarks of species q in resonance i, and the chemical 
potentials correspond to those in the QGP. We take an excluded volume, v exc = 4nRl xc /3, 
for every resonance in the RG phase ||. Our justification for this is that when the resonances 
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overlap, the region where this occurs should be in the QGP phase, so we do not allow overlap 
in the RG phase. 

The transition temperature, T c , is obtained by setting the pressures equal in the two 
phases, with all chemical potentials fixed. In Fig. 1, we show the transition temperature as 
a function of fi B , for the cases (i) n s = rig, fi-g = —fig (strangeness equilibrium under the 
strong interactions), and (ii) fi-g — fis — > — 00 (complete strangeness suppression). We give 
results for R exc = 0.5 and 1 fm in each case, adjusting the bag constant to give T c = 150 MeV 
for fi B = 0; for the equilibrium case, we take 5 1 / 4 = 217 and 220 MeV for R exc = 0.5 and 1 
fm respectively, while for the case of complete strangeness suppression we take B 1 ^ 4 = 207 
and 210 MeV. Here (and for the remainder of this paper) we evaluate all momentum integrals 
numerically to one percent accuracy. 

We fix /is for the strangeness equilibrium curves by requiring that the QGP or RG be 
strangeness-neutral (rij = n s ). At fixed fi B , the transition from one strangeness-neutral 
phase to the other phase (of arbitrary strangeness) occurs at almost exactly the same tem- 
perature whether the initial phase is QGP or RG, so we only show RG curves. This result, 
first noted in Refs. PJIQ"|, is somewhat surprising, as the value of fis in the RG depends 
on the strange hadron spectrum, so that the transition temperatures could easily be very 
different for the two phases. 



III. EVOLUTION OF THE HOT MATTER 



We model the initial evolution following Ref. [14]]. We approximate the nucleon wave- 
functions to be constant inside cylinders with radius tn and length (parallel to the beam 
direction) In, and zero outside these cylinders. The radius is given by the nuclear radius, 
tn = 7 fm for Au, while the length in the center of momentum (cm) frame is In ~ ?"jv/7- 
Here 7 = (1 — w 2 ) -1 / 2 , where v « 1 is the nuclear velocity in the cm frame; for fixed target 
collisions at beam energy 11.6 GeV/nucleon, 7 = 2.7. In principle, In is bounded from below 
by the minimum parton wavelength, of order m^ 1 = 1.4 fm, but for the collisions considered 
here In = 2.6 fm so this lower bound is unimportant. 

We then assume that the nucleon properties are unchanged during the collision, and 
neglect collisions of secondary particles. The rate of production of any quantity during the 
collision is then proportional to the overlap of the nucleon wavef unctions. For example, if 
the cylinders first touch at proper time r = \/t 2 — z 2 = 0, where z is the position along the 
beam axis, the baryon rapidity density at rapidity y = is 

f 0, r* < 0, 

dN B ,dy = k B l% x f_ 2(1 _ n)2] A J < £ < VJ, (2) 
{ 1/2, 1 < r,, 

where k B is an unspecified normalization constant, and t* = r/Zjy- The volume of hot matter 
per unit rapidity is approximately vrr^r, so the baryon density at y = is 
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0, r* < 0, 

dN B /dy I 2r*, < r* < 1/2, 

^ Trr^Zjv X | [1-2(1 -T*) 2 ]/r*, 1/2 < r* < 1, 

1/r*, 1 < r*, 



where dN B /dy is the observed baryon rapidity density at ?/ = 0. 

For r > /^v, we evolve the hot matter hydrodynamically, maintaining boost-invariance 
(and hence neglecting transverse expansion). We assume strong superheating in the RG 
phase, so that the hot matter never makes a transition to QGP, as the temperature never 
rises very far above T c . The hot matter then follows the boost-invariant hydrodynamic 



equations [15 



de -(e + P) 



dr t 
driB — tib 
dr t ' 
d(n s - n?) 



(5) 
(6) 



dr 

along with the equilibrium condition, 

= -t*s- (7) 

Here e, P, and rig are respectively the energy density, pressure, and s density. We connect 
the hydrodynamic evolution to the initial evolution by using the fact that the baryon and 
entropy rapidity densities, respectively dNB/dy and dS/dy, are approximately invariant after 
the initial collisions. 

In Fig. 2, we show the trajectory followed by the hot matter at y = in a central 
Au+Au collision at 11.6 GeV/nucleon, using Rf = 0.5 fm. We use the preliminary results 
dNs/dy = 120 and dN n + +n - jdy = 120 [|T^] in the central rapidity region, from which we 
estimate dS/dy ~ 1700 if the freezeout temperature Tf = 100 — 120 MeV. At r = In/2 
and t = In, T = 162 MeV, and the maximum temperature reached is 169 MeV, so the hot 
matter is only superheated by 20 — 30 MeV during the initial evolution. For comparison, we 



show also the trajectory obtained with the initial conditions used in Ref. [IT|, e ~ 2 GeV- 
fm~ 3 and % ~ 0.8 fm -3 for r < In, along with the critical temperature, T c . The maximum 
temperature reached with this second trajectory is 200 MeV, significantly higher than for 
our trajectory. However, our trajectory should reflect the conditions of the thermalized 
matter, which is most likely to remain at y — 0, and thus will be accurately determined by 
rapidity density measurements. 



IV. DECAY TO OBSERVED LEPTON AND KAON PAIRS 

There are a number of calculations of the <fi mass and width at finite T. The range of 
predictions for the mass is very wide. Haglin and Gale find that the mass shift Sm(T) 
increases slowly and monotonically with increasing T, with <5m(200 MeV) ~ 4 MeV, while 
Asakawa and Ko find that 5m decreases monotonically with 5m (190 MeV) ~ —170 
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MeV. The range of values for the width, T(T), is much narrower. Ko and Seibert Q and 
Haglin both find that the RG contribution to T(T) is approximately 30 MeV at T = 200 
MeV, although it may be as low as 10 MeV if vertex form factors are included H. The QGP 
contribution to the width depends on the dynamics of the two-phase coexistence region, as 
it is proportional to the fraction of matter in the QGP phase and inversely proportional to 
the mean QGP droplet radius ||. In the adiabatic limit (when the nucleation rates for both 
phases are infinite), the QGP contribution to T(T) is infinite; however, the contribution to 
observable quantities is probably small, just a few MeV in an Au+Au collision at yfs = 200 
GeV/nucleon and of order 10 MeV at ^Js = 10 — 20 GeV/nucleon. 

For the mass, we take two very different parametrizations, following Refs. H and H; 



m 



(T) = 1020 + 4(T/200 MeV) 4 MeV, (8) 
m^ K {T) = 1020 - 200(T/200 MeV) 4 MeV. (9) 

We use these masses in the RG phase, and assume that there are no recognizable <j) mesons 
in the QGP. As the corrections for finite baryon density are uncertain, we neglect them at 
present. For the width, we use 

T(T) = 4.43 + 25(T/200 MeV) 2 MeV. (10) 

The dilepton signal from before freezeout is obtained by convoluting the four-volume per 
unit rapidity, 7rr 2 Tdr, the </> density n^T), the thermal phi mass distribution p(m, T), and 
the decay rate to dileptons, IV The post-freezeout contribution is calculated in the same 
manner, except that we use the vacuum mass and width, and assume that after freezeout 
the number of mesons decreases exponentially in r with the T = decay rate. 



dm 2 dy 



T f n^(T f )p(m,0) 



r(o) Jt +i n /2 



+ / dr t n$(T) p(m, T) 



(11) 



We take Tu to be constant, evaluate numerically, and assume a Breit-Wigner mass dis- 
tribution, 

P(m, T) = m ^ r(T)/7r . (12) 

m 2 - mJ(T)l + m 2 (T) T 2 (T) 

We neglect possible dileptons produced before the hot matter is thermalized at proper time 
ro ~ 1 — 2 fm/c after the point of maximum overlap; this signal is probably small, and we 
have no idea what the mass distribution is for these particles. The calculation is simplified 
by the fact that the lepton mean free path is long, so that virtually all produced pairs escape 
from the hot matter without interacting. 

The kaon signal is calculated in almost the same manner, but with two changes. First, 
the decay rate to kaons, Tkk, depends strongly on the mass. We include this by using 

r„cn - ^^r(o), (13) 

m^(T) 2 K(0)2-4m|.] 3/2 
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where tuk = 495 MeV is the vacuum kaon mass. Second, we fold in the probability for the 
kaons to escape the hot matter without interaction, which is approximately 

= 2,r N r x r K (T)v4T)/S 
T K (T)vt(T) 

where tk is the mean time between kaon interactions, and v$ is the mean thermal velocity 
for cf) mesons. The numerator in the first expression is the product of the surface area of 
the cylinder and the mean distance from the surface for a pair of emerging kaons (half that 
for a single kaon), while the denominator is the total volume. As the kaon velocity in the 
rest frame is not too large, we assume that the kaon and <j) velocities are the same; when 
the velocity is small enough that this is not true, it is very unlikely that both produced 
kaons will escape without interaction, so this should produce only a small correction since 
our escape probability is proportional to and hence goes to zero in this case. 



dN KK 2 



dm 2 dy 



nr 



Tf7if(Tf)p(m,0) + / dTTn^(T)T KK (T)VKK(T)p(m,T) 

Jto+In/2 



We take to be the (non-relativistic) rms thermal velocity, 



and use the parametrization 



3T 

m <t> (T)_ 



1/2 



• (16) 



(17) 
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/ lSOMeV V , 
r K = 3 ^ J fm/c, (18) 



which agrees well with the results of Haglin and Pratt JT7 



We show predicted dilepton and dikaon mass spectra in Fig. 3, taking r = and Tf = 100 
MeV. The dikaon spectra differ significantly only far from the peak, when the signal has 
dropped to less than 1% of the peak value. As there is large combinatoric background near 
the peak, the small difference in the dikaon signal will almost certainly be lost, while the 
difference in the dilepton signal is over 10% and thus may be detectable. 

Note the small secondary peak in the dilepton spectrum at m ~ 920 MeV for m§ K 



This peak appears because the hot matter spends 1 — 2 fm/ with T ~ 160 MeV, just as in 
the case of a first-order phase transition |2||| (when the hot matter remains at T c for a long 
time). However, the peak here reflects the (approximately constant) initial temperature, 
which results from the spatial extent of the parton wavefunctions, and depends only weakly 
on T c . In the case of a phase transition, the temperature of the hot matter does not remain 
constant in the coexistence region even though T c is the same for the two phases, because 
the entropy per baryon is significantly higher in the QGP at fixed T and |19| |, leading to 



smearing of any secondary (p meson peak created during a phase transition at finite //g. The 
larger peak at m ~ 1000 MeV is due to the fact that the system is assumed to freeze out 
instantaneously, so that there is no hot matter with < T < Tf and thus no <fi mesons are 
seen in the corresponding mass range; this peak will soften and possibly disappear if more 
realistic dynamics are used. 
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V. CONCLUSIONS 



Although the total signal from dikaons is larger than that from dileptons, the vast ma- 
jority of the dikaons come from mesons that decay after freezeout of the hot matter. Thus, 
we see that, not very surprisingly, it is difficult to measure thermal mass shifts by observing 
strongly interacting particles. Conversely, dileptons make a good probe of thermal mass 
shifts, since they escape from the hot matter throughout the collision. 

There are two reasons for the lack of a strong signal from dikaons. The first is that, 
once the mass drops appreciably, the decay rate to dikaons vanishes. However, this is not 
so easily corrected by observing meson decay to lighter particles, such as pions. Although 
the decay rate then increases, the velocity of the decay products also increases, so that it is 
greater than the mean thermal velocity of the heavy meson. The decay products then tend 
to emerge back-to-back in the lab frame, so that the probability of both decay products 
escaping from the hot matter without interacting is small. It is also possible that the kaon 
mass drops with increasing T |[20| , and that the mesons decay to these lighter kaons. 
However, in this case we expect that the thermal kaons will interact when they leave the hot 
matter, as they acquire the vacuum kaon mass, and that this interaction will likely remove 
the pairs from the peak. 

As a result of the lack of sensitivity of the dikaon signal to thermal meson properties, 
it is very important to study dilepton production in ultra-relativistic heavy ion collisions, 
particularly near vector meson peaks. This is the only likely means for detecting finite- 
temperature shifts in meson masses and widths. As the study of meson properties may pro- 



vide valuable insight into the equation of state of hot strongly-interacting matter |21|J2^ , 
we would hope to see more experimental attention devoted to these measurements in the 
future. 
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FIGURE CAPTIONS 



Fig. 1: T c vs. hb using R exc = 0.5 and 1 fm, for the cases of strangeness equilibrium under 
the strong interactions and of complete strangeness suppression (denoted "no s"). 

Fig. 2: Trajectories followed by the hot matter produced in a central Au+Au collision at 
lab energy 11.6 GeV/nucleon. The solid trajectory follows the evolution described in 



the text, while the KVV trajectory uses the initial conditions from Ref. |LT| as the 
starting point for the hydrodynamic evolution. The matter is assumed to remain in 
the RG phase, with R eX c = 0.5 fm. For comparison, T c is shown for a chemically 
equilibrated resonance gas with n s = rig, taking B = (217 MeV) 4 . 

Fig. 3: Predicted dikaon (kk) and dilepton (11) spectra from a central Au+Au collision at 
lab energy 11.6 GeV/nucleon, using m^(T) from Refs. and 0, denoted by HG and 
AK respectively. Evolution of the hot matter proceeds as described in the text, with 
Rexc = 0.5 fm. 
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Fig. 1 - Seibert 
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Fig. 2 - Seibert 
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